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DIFFUSION  FIELDS  IN  A  TOROIDAL  CONDUCTING  SHELL 
OF  CIRCULAR  CROSS  SECTION 


I.  Introduction 

Over  the  years  abundant  analytical  work  has  been  done  on  the  diffusion  of  the  mag¬ 
netic  field  in  hollow  circular  cylindrical  conductors  of  infinite  length1-3,  and  has  led  to 
the  understanding  of  the  diffusion  process  in  such  a  geometry.  Much  less  attention  has 
been  given  to  the  diffusion  of  fields  in  toroidal  vessels.  As  a  rule,  the  diffusion  problem 
through  toroidal  chambers  is  more  complex  than  through  cylindrical  conductors.  The  dif¬ 
ficulty  in  toroidal  geometry  is  due  mainly  to  the  application  of  the  boundary  conditions 
in  the  inner  and  outer  surfaces  of  the  torus.  However,  when  the  wall  thickness  of  the 
torus  becomes  small  in  comparison  to  its  minor  radius,  i.e.,  for  a  conducting  shell,  then 
the  problem  simplifies  considerably.  The  toroidal  shell  approximation  has  been  used  to 
compute  analytically  the  magnetic  field  of  two  toroidal  segments,  that  carry  known  time 
independent  currents.  In  this  paper,  we  have  extended  the  work  of  P.  Rolicz  et  al.4  to 
study  the  diffusion  of  fields  through  toroidal  conducting  shells. 

The  analytic  results  are  further  simplified  under  the  additional  assumption  of  a  small 
aspect  ratio  toroidal  conducting  shell.  To  the  lowest  order  on  the  aspect  ratio,  our  work 
indicates  that  the  diffusion  of  fields  through  a  conducting  shell  of  conductivity  cr,  major 
radius  r0,  minor  radius  a  and  thickness  dw,  depend  on  three  characteristic  times  t0,  tx 
and  r2,  where  tx  =  n0oadw/ 2  (diffusion  time),  r0  =  2rj[£n(8r0/a)  -  2]  (L/R  time),  and 
r2  =  t'i/2.  In  contrast,  the  diffusion  process  in  cylindrical  shells  depends  on  a  single 
characteristic  time  (rj). 

It  is  well  known  that  there  is  a  time  lag  between  the  magnetic  field  in  the  interior  of 
the  conducting  shell  and  the  externally  applied  magnetic  field.  This  time  lag  can  be  easily 
measured  in  the  laboratory  and  therefore  is  a  useful  physical  quantity.  When  the  transient 
phase  of  the  diffusion  process  is  over,  the  time  lag  becomes  independent  of  time  and  is 
called  the  delay  time  .  For  a  linearly  rising  applied  magnetic  field  the  delay  time  for  a 
cylindrical  conducting  shell  is  equal  to  T\.  On  the  other  hand,  in  a  toroidal  conducting 
shell  the  delay  time  depends  on  both  To  and  T\  and  also  on  the  flux  condition,  i.e.,  the 
ratio  of  the  average  field  to  the  local  field  on  the  minor  axis.  Specifically  =  (r0  +  3rx)/2 
for  the  special  case  of  a  field  that  satisfies  the  betatron  flux  condition,  i.e.,  when  the  value 
of  the  flux  condition  is  equal  to  two. 

In  contrast  to  the  field,  the  current  flowing  on  the  wall  of  the  conducting  shell  depends 
only  on  to.  Thus,  by  measuring  experimentally  the  delay  time  of  the  fields  and  also  the 
rise  time  of  the  wall  current,  all  characteristic  times  can  be  determined. 
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The  results  of  the  analysis  have  been  compared  with  the  predictions  of  the  two- 
dimensional  finite  difference  numerical  code  TRIDIF5.  The  agreement  between  theory 
and  numerical  computation  is  fairly  good,  in  spite  of  the  fact  that  the  theory  was  done 
for  a  toroidal  conducting  shell,  while  in  the  numerical  computation  the  torus  had  finite 
thickness.  Specifically,  the  delay  time  predicted  by  the  analysis  is  approximately  5%  longer 
than  the  prediction  of  the  code.  Similarly,  the  temporal  profile  of  the  wall  current  from 
the  analysis  and  the  TRIDIF  code  agree  to  within  a  few  percent.  However,  the  temporal 
profile  of  the  external  field  index  and  flux  predicted  by  the  theory  and  the  code  are  in  good 
agreement  after  a  few  delay  times  but  not  initially. 

The  delay  time  predicted  by  the  theory  and  the  code  has  been  compared  with  the  delay 
time  measured  in  the  toroidal  chamber  of  the  NRL  modified  betatron.  For  the  parameters 
of  the  NRL  device,  the  computed  delay  time  is  37  /zsec,  while  the  measured  delay  time  is 
34  /usee.  The  less  than  10%  difference  between  the  two  delay  times  is  probably  within  the 
uncertainty  of  the  measurement. 

In  Section  n,  the  theory  is  developed  and  a  comparison  is  made  with  the  numerical 
results  obtained  from  the  TRIDIF  code.  Section  III  contains  a  description  of  the  NRL 
toroidal  device  and  of  the  coils  that  generate  the  external  magnetic  field,  and  the  measured 
quantities  are  compared  with  the  theoretical  results.  Finally,  in  Section  IV,  the  conclusions 
are  presented. 


II.  Theory  and  Comparison  with  Numerical  Results 

The  diffusion  fields  are  computed  near  the  minor  axis  of  the  circular  cross  section, 
toroidal  conducting  shell.  As  shown  in  Figure  1,  the  conductor  has  a  major  radius  r o,  a 
minor  radius  a,  wall  thickness  dw,  where  dw  «  a,  and  conductivity  a.  In  the  presence 
of  an  external  magnetic  field  which  is  axisymmetric  and  time-dependent,  the  diffused 
magnetic  field  inside  the  conductor  is  described  by  the  diffusion  equation 

dA 

VxVxA  =  -noo—,  (1) 


where  the  vector  potential  A  has  only  one  nonzero  component  Ap  that  depends  only  on  the 
cylindrical  components  (r,  z)  and  on  the  time.  The  magnetic  field  components  are  given 
by 


Br  =  - 


dAp 

dz 


„  1  drAp 

r  dr 


(2a) 

(26) 
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while  the  electric  field  is  equal  to 


Ee  =  -• 


Equation  (1)  is  identical  to  Ampere’s  law  combined  with  Ohm’s  law  inside  the  con¬ 
ductor.  In  the  special  case,  when  the  toroidal  conductor  is  a  shell,  i.e.,  when  dw  <<  a, 
this  equation  can  be  integrated  and  provides  the  first  boundary  condition,  i.e., 

n  x  (Bout  -  Bin)  =  n0odwEin ,  (4) 

where  n  is  the  unit  vector  normal  to  the  conducting  wall  and  directed  towards  the  region 
outside  the  conductor  and  no  is  the  permeability  of  the  vacuum.  The  second  boundary 
condition  is  obtained  from  the  requirement  that  the  electric  field  is  continuous  across  the 
boundary,  i.e., 

Ein  =  Eout.  (5) 

For  a  small  aspect  ratio  a/ro  vessel,  it  is  appropriate  to  use  toroidal  coordinates.  The 
cylindrical  coordinates  (r,  z)  are  related  to  the  toroidal  coordinates  (77,  £)  by: 


r  =  b 


sink  rj 

cosh  77  —  cos  £  ’ 


z  =  b  v 

cosh  77  —  cos  £  ’ 

where  b  is  a  constant.  These  relations  can  be  easily  inverted,  namely: 


■2„_  (riY  _  (r-b)*+z* 
\^2  )  (r  +  6)2  +  z 2  ’ 


e  n  cos 


From  Eq.  (7a),  we  see  that  for  fixed  77,  the  coordinates  (r,  2)  describe  a  ::  ;cle  whose  radius 
is  bj  sinhT7.  If  for  77  =  770  this  circle  coincides  with  the  toroidal  shell  whose  minor  radius  is 


a,  then  it  is  straightforward  to  show  that  6  =  roj^l—  ^  —  JJ  .  The  points  (77,  £)  outside 

the  toroidal  shell  are  determined  by  the  inequality  77  <  770,  while  the  points  inside  the 
hollow  region  of  the  shell  are  determined  by  77  >  770.  In  both  regions,  the  right  hand  side 
of  Eq.  (1)  is  zero  and,  in  toroidal  coordinates,  this  equation  reduces  to 


-te)'] 


hrjh^he  [dr) 


C^jrir) +  rf  {hjr-d-jf)  1  -  hA‘  =  °- 

V  hr,  drj  )  di  \  h£  d£  )  J  h2 
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where 


(9a) 


hr, 


he  = 


-  6 

cosh  77  —  cos  £  ’ 

6  sinh  77 

cosh  rj  —  cos  £  ’ 


(96) 


If  we  set  Ae  =  6(cosh  —  cos  £)  *  F(tj,  £),  then  the  differential  equation  satisfied  by  F(tj,  £)  is 
a  separable  equation  and  its  solution  can  be  expressed  in  terms  of  the  toroidal  functions. 
Specifically,  the  solutions  of  Eq.  (8)  in  the  two  regions  inside  and  outside  the  toroidal  shell 
are: 

OO 

Agn  =  6(cosh  rj  -cos£)*  ^  emam (cosh  77)  cos  m£,  (10a) 

m= 0 


00 

A%ut  =  Aeext  +  6(coshrj  —  cos£)s  ^  em6m(t)Pt^_i(cosh»7)cosm^,  (106) 

m=0 

where  eo  =  1,  £m  —  2  when  m  —  1,2,3,...,  P^-i  (cosh  r]),Qlrn_L  (cosh  77)  are  the  associated 
Legendre  functions  of  the  first  and  second  kind,  respectively,  and  Agxt(r, 7,  £,  t)  is  the  vector 
potential  associated  with  the  applied  external  field.  The  time-dependent  parameters  am(f) 
and  6m(t)  are  determined  from  the  boundary  conditions.  Since  the  ^-component  of  the 
magnetic  field  is  equal  to 

Bi  =  ~h^hed^heAe' 

and  the  external  field  is  zero  at  t  =  0,  the  boundary  conditions  of  Eqs.  (4)  and  (5), 
expressed  in  terms  of  thp  vector  potential  A$,  become: 


A'en(rjo,(,t)  =  A%ut(r?0,{,t), 


(12  a) 


,  ,  /AF'in.C'tl-Arh'C't) 

(cosh  770  -  cos  0 2  5-  (  — r - ~rri - 

or\  \  (cosh  77  —  cos  f)  2 


) 


—  T 


dA?{  no^,t) 

dt 


(126) 


77=770 


where  r  =  fi0odwb. 

In  the  following,  we  shall  assume  that  there  are  no  external  current  coils  very  near 
the  toroidal  conductor  or  inside  it.  If  r]c  is  the  toroidal  coordinate  of  the  coil  nearest  to 
the  minor  axis,  and  such  that  t]e  <  770,  then  the  external  vector  potential  Agxt  is  given  by 


OO 

Aeext  =  6(cosh  77  -cos£)’  JZ  €m«m*(0Gm_$  (cosh?)  cos  m£,  (13) 

m= 0 

for  all  77  >  77c  The  time-dependent  parameters  a”*(i)  are  known  and  are  associated  with 
physical  properties  of  the  external  field,  such  as  external  field  index,  flux  condition,  etc. 
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Now,  it  is  straightforward  to  compute  the  unknown  coefficients  am,bm  from  the  bound¬ 
ary  conditions.  From  Eq.  (12a)  we  obtain 


b  m  = 


{dm  -  (cosh  170) 

■^m— J-(C°sh  T?0) 


Similarly,  Eq.  (12b)  together  with  Eq.  (14)  and  the  identity6 

'  du  }  du  4(u2  —  1)  ’  (  ’ 

leads  to  the  following  set  of  coupled  differential  equations  for  om: 

1  1  1  Ql(coshr?0) 

do  = - (00  -  a‘0xt)  +  —  r — — Z — \  (a*  “  •“*) .  (16a) 

To  t0  cosh  rjo  (cosh tj0) 


1  /  e  t\  1  1  Qlm-  |(COshr?o) 

am  = - {am  -  a^)  + - - - - — — J — — - - - r  (am_!  -  ami1; 

Tm  rm_!  2coshr7o  Q^_i(cosh»7o) 


,_1 _ 1  ei,+i  (cos^o)  «  . 

rm+1  2coshi)0  <?’ _j.(cosh>)o)  m+'  m+1  ’ 


(16*) 


where  m  =  1, 2, 3, ...  and 


^m2- §  (c°sh  1  (cosh  *»o)  • 


The  toroidal  functions  and  appearing  in  Eqs.  (16a),  (16b)  and  (17)  are 

given  by  the  following  exact  expressions6: 

P"-t(cosh.)  =  -t- 


(-l)‘2M1r(m  +  »tl)  lmWil,  y-  (n+l),(m  +  i.+  l), 

7rlr(m-|-l)  JZ'q  s!(m  +  l). 


*  [£n(4en)  +  Ua  +  Um  +  8  Vn-(-a  ^m+n+«  > 


(18a) 


5 


(186} 


i)(sinh,)ne— (— i)n 


y-  {n  +  ^)a{m  +  n  +  ^)e  _2arj 

k  *!(m+i)a  e  ’ 


where 


^  s!(m  +  l)a 

(a)fl  =  a(a  +  l)(a  +  2). ..(a  +  s  —  1),  (°)o  =  1» 

1  n  1 

Un=2^k'  u°  =  0’ 

Jt=i 

n  1 

’'»  =  Esr^T'  Vo=0- 


(19a) 

(196) 

(19c) 


For  m  —  0,  the  first  term  in  Eq.  (18a)  is  omitted,  since  6mn  =  1  for  m  =  n,  and  6mn  =  0  for 
m  7^  n.  The  expressions  above  are  appropriate  for  the  region  inside  the  toroidal  conductor, 
i.e.,  for  T)  >  770,  as  well  as  on  its  surface. 

Up  to  this  point  the  results  are  exact.  For  small  aspect  ratio  a/ro,e~n°  =  (ro  +  a  — 
6)/(r0  -h  a  +  6)  «  a/2r0  <<  1.  Thus,  to  the  lowest  order  in  the  aspect  ratio,  Eq.  (17) 
becomes 

T0  =  2td  £rt“  -  2  ,  (20 a) 

a 


'  m  —  ) 

m 


where  m  =  1,2,3, ...,  and 


td  = 


/^0  @d\jjCL 


(206) 

(21) 


The  time  To  is  the  L/R  time  of  the  toroidal  conductor,  while  rm  are  the  higher  order 
diffusion  times. 

An  approximate  expression  for  A'6n  can  be  obtained  from  Eq.  (10a)  by  keeping  terms 
up  to  order  e~2rl,  i.e.,  up  to  second  order  toroidal  corrections.  Making  use  of  the  identity6 


(cosh  r?  —  cos  £)  *  =  £mCm(t))  cos  m£, 


where 


(r))  cos  m£, 

(22) 

le-*  1 

1  +  e-2^(coshr?)j’ 

(23a) 

2”) 

(236) 

6 


*  Qm-  i  (cosh  r?)  -  1  *e-2n  (fim- 1  (cosh  7?)  +  Qm+ 1  (cosh  r})) 
and  also  of  Eq.  (18b),  a  lengthy  calculation  leads  to  the  approximate  expression 


7T 


aQ+-{ao-ai)e  2" 


—  (ao  —  3ai)e  n  cos  £  —  -{a0  +  6ai  —  15a2)e  2r,cos2£ 


(24) 


Using  Eqs.  (7a),  (7b)  and  keeping  terms  up  to  order  ^1  —  and(z/r)2 ,  Eq.  (24)  reduces 

(25) 


to 


4*n 

Ae 


2V2 


°o  ~  ^(a0  -  3ai)(l  -  £) 


1  b  3  z 

+  — (a°  +  15a2)(l - )2  +  —  (ao  +  4ai  —  5a2)(-)2 

lb  r  lb  r 


i.e.,  the  approximate  A™  has  been  expressed  in  terms  of  the  cylindrical  coordinates  (r,  2). 
The  magnetic  field  components  near  the  minor  axis  are  determined  from  Eqs.  (2a)  and 
(2b)  and  are  equal  to 

(26a) 


B\n  «  - Bzono  -, 


B\n 


Bgo 


where 


Bzo  = 


1  -  Up(  1 - ) 


(ao  +  3aj), 


4v/2 


is  the  vertical  field  component  near  the  minor  axis  and 


3  ao  ■+■  4ax  —  5a2 

nD  = 


(26  b) 
(27a) 

.  ,  .  (27 b) 

4  ao  +  3ai 

is  the  field  index.  To  first  order  in  a/ro,  the  magnetic  flux  $0  through  the  area  7rrQ  on  the 
midplane  is 

$0  =  2nr0A'9n(r0,0,t)  (28) 

“  -2’rl^ao 

Defining  ep  by  means  of  the  relation  $0  =  2irrQBzo(l  —  tp),  the  vertical  component  of  the 
field  on  the  minor  axis  can  be  written 


Bzo(l  —  ep) 


7T 


2>/2 


ao. 


(29) 
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The  parameter  cr>  is  a  measure  of  the  flux  condition.  When  to  =  0,  the  betatron  flux 
condition  is  satisfied. 

The  quantities  ao,Gi,a2  are  determined  completely  by  the  magnetic  field,  the  field 
index  and  the  flux  on  the  minor  axis,  and,  therefore,  are  parametric  representations  of 
these  physical  quantities  But  the  vector  potential  that  describes  the  external  field  in  Eq. 
(13)  has  exactly  the  same  form  as  Axen  in  Eq.  (10a).  Therefore,  if  Bo,  n  and  $QXt  = 
27rroi?o(l  —  f)  are  the  external  field,  external  field  index  and  external  flux  on  the  minor 
axis,  then  a,QXt,  a\zt  and  a\xt  can  be  expressed  as  follows: 


7T 

(30a) 

°r'  =  “17Bo(1  +  £)' 

(306) 

1 D7T 

(30c) 

These  expressions  are  used  in  solving  Eqs.  (16a),  (16b).  It  is  assumed  that  only  B0 
is  time-dependent.  This  is  a  valid  assumption  provided  that  the  coils  that  generate  the 
external  field,  carry  the  same  current  and  are  filaments,  i.e.,  there  is  no  diffusion  process 
associated  with  them. 

All  the  physical  quantities  of  interest  have  been  computed  to  first  order  in  toroidal 
corrections.  Therefore,  it  is  sufficient  to  compute  the  parameters  ao,ai,a2  to  first  order 
in  a/r0.  If  all  terms  of  order  e~2r>0  or  higher  are  omitted,  Eqs.  (16a),  (16b)  become 

d0  «  -  —  (oo  -  ae0xt),  (31a) 

To 


~~{am  —  am)  + 


2m 


(am-  1  am— 1)> 


(316) 


Tm  Tm— j  2m  4~  1 

where  m  =  1,2.  This  system  of  equations  can  be  integrated  easily.  Under  the  initial 
conditions  am(0)  =  0,  the  solutions  for  r2  =  T\j2  are: 


o0 


2\/2 

7T 


(1  -  e)A0, 


(32a) 


a  i  = 


a2  =  - 


2v/2 
3x 

2y/2 

15jt 


(' +  e  -  r^(1 "  £>) Ai  +  tH{x  ~ e)Xo]  ’ 

(  7  —  8n  +  e  —  4(l-fe)  +  — - (1  —  e;  J  A2 

L  V  2  - a  / 


(32  b) 
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+  (“t1  +  £)  -  rr^1  -  £>)  +  (1  —  a) (2  —  a) (1  -  £)Xo.  1 


where 


Am  =  — e-*/T"  f  (336) 

JO 

Substitution  of  ao,ai,a2  into  Eqs.  (27a),  (27b)  leads  to 

B'° =  (5 +  rh;) (1  - e)Ao  + 11  “  G  +  rrs) (1  -  £)] Au  (34o) 

(§  +  *)(»- e Mo  +  ["-(§  +  A)d- £)] A* 

nD  =  7 - r - - 7- - \ -  •  (346) 

(5  +  if;)  (1  -  £Mo  +  [»  (5  +  A)  (!  "  ')]  A> 

Another  interesting  physical  quantity  is  the  surface  wall  current  density,  which  is  given 


(32c) 


(33a) 


(336) 


(34a) 


(346) 


J ad  —  od\ 


dAr{fio,t,t) 


To  first  order  in  toroidal  corrections,  Eqs.  (25),  (31a),  (31b),  (32a),  (32b)  give 


J8e  «  ^-{a(l  -  c)(A0  -  Bo)  -  (4  +  a-- - (1  -  e)(A0  -  B0)  (36) 

fjL0a  ^  L\  I  -  a/ 

-  +  6  ~  “ e))  G*1  _  5°)  ^cos^}’ 

where  cos  4>  =  (r  -  r0)/a  and  r  is  the  radial  position  of  any  point  on  the  wall  chamber. 
The  wall  current  is  equal  to 

1=  [  Jsead(f>  (37) 

Jo 

=  —a(l  -e)(A0-  Bo). 

Ho 

As  a  first  application,  consider  a  rising  external  magnetic  field  that  reaches  a  constant 
value,  i.e., 

B0(£)  =  60(1  -  e_t/r<!1*).  (38) 


From  Eq.  (33b),  we  have: 


__  L  1  1  -  — t/ Text  ,  TmlUxt. -t/rm 

~b°  1  ~  1  rm  e 

Text  TeIt 
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For  t  >>  Text  and  t  >>  rm,  the  coefficient  Am  «  bo  and  from  Eqs.  (34a),  (34b),  Bzq  «  b0 
and  np  w  n.  Also  the  wall  current  goes  to  zero.  All  three  conclusions  are  as  expected. 
For  a  sinusoidally  varying  field 


Bo(t)  =  bo  smut. 


(40) 


the  coefficients  Am  are  given  by 

bo 


Am  — 


1  +  (wrm)2  L 


smut  —  urrn  cos ut  +  uTme  t^T" 


(41) 


After  a  long  time,  i.e.,  for  t  >>  rm,  the  external  field  and  the  field  on  the  minor  axis  are 
not  in  phase,  but  their  phase  difference  remains  constant.  On  the  other  hand,  for  small 
times,  i.e.,  for  ut  <<  1  and  t  «  rm,  the  coefficients  Am  «  b0ut2 /2rm  and  the  field  index 
becomes 


np(t  =  0) 


16n  —  (5a  +  6)(l  -  e) 
4[2-(l  +  a)(l  — e)l  ' 


(42) 


Furthermore,  assuming  that  urm  <<  1  and  ut  <<  1,  but  not  t  >>  rm,  we  obtain 


A 


m 


b0u  -  rm(l  -  e  , 


(43) 


and 


Bzo  «  bo 


t-Td  + 


T-*— 

To  —  Ti  TO  ~  T1 


where  the  delay  time  ra  is  equal  to 


Td  = 


(1  -  e)r0  +  (3  -  e)ri 
2 


(44) 


(45) 


In  v  ntrast,  the  delay  time  for  a  cylinder  of  infinite  length  is  equal  to  the  diffusion  time. 
It  is  apparent  from  Eq.  (45)  that  in  the  case  of  the  toroidal  conducting  shell,  the  delay 
time  depends  linearly  on  both  the  L/R  time  r0  and  the  diffusion  time  T\.  In  addition,  for 
a  toroidal  shell,  Bz o  has  an  exponentially  decreasing  dependence  on  both  times  To  and  T\, 
while  the  wall  current,  from  Eq.  (37),  depends  only  upon  r0,  i.e., 

/  =  1  -  c)  (\  -  ,  (46) 

Mo  v  ' 


i.e.,  it  reaches  exponentially  its  maximum  value  with  the  L/R  time  To. 

A  direct  comparison  has  been  made  between  the  theory  presented  above  and  the 
numerical  results  from  the  2-dimensional  TRIDIF  code.  The  set  of  parameters  used  for 
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the  comparison  axe  listed  in  Table  I.  For  these  parameters,  Eqs.  (20a),  (20b),  (45),  give 
r0  =  11.28  n sec,  T\  =  2.95  /zsec,  a  =  0.26,  and  r &  =  10.42  /*sec.  The  solid  curves  in  Figs. 
2(a)-2(d)  show  Bz o,  np,  $o/{’jrroBzo),  and  7  as  a  function  of  time  .  In  the  computer  run, 
the  inner  minor  radius  of  the  toroidal  conductor  is  15  cm  and  the  outer  minor  radius  is  17 
cm.  In  the  region  of  the  chamber,  the  mesh  has  a  cell  size  A r  =  0.5  cm,  A z  =  1.0  cm,  but 
further  out  it  is  much  less  dense.  Thus,  in  the  radial  direction  there  are  157  mesh  points 
over  a  distance  of  300  cm,  while  in  the  vertical  direction  there  are  71  mesh  points  over  a 
distance  of  200  cm,  and  the  time  step  used  is  A t  =  0.5  jzsec.  Advantage  is  taken  of  the 
midplane  symmetry,  i.e.,  the  calculation  is  confined  to  the  region  z  >  0.  The  external  field 
is  generated  by  current  filaments  (coils)  located  at  the  same  positions  as  in  the  modified 
betatron  experiment.  All  the  coils  carry  the  same  current  and  have  the  same  sinusoidal 
dependence.  In  the  absence  of  the  toroidal  chamber  the  maximum  vertical  field  generated 
by  these  coils  on  the  minor  axis;  i.e.,  at  r  =  100  cm,  z  —  0,  is  Bo  =  1267  Gauss,  when  the 
maximum  current  through  the  coils  is  40  kA.  This  is  the  value  of  60  used  in  the  theory.  In 
the  absence  of  the  toroidal  chamber,  these  coils  generate  a  field  index  of  0.54  on  the  minor 
axis  and  an  e  =  —0.067,  but  the  values  chosen  in  the  theory  are  slightly  different,  so  that 
at  times  t  >  30  fisec  the  theoretical  and  numerical  nD  and  $o/(nroBzo)  are  identical. 

The  dashed  curves  in  Figures  2(a)  -  2(d)  give  the  predictions  of  the  code  for  Bz0,  np, 
$o/{nrlBzo)  and  7  as  a  function  of  time.  The  straight  line  in  Fig.  2(a)  is  Bz0,  when 
a  =  0.  The  numerical  delay  time,  obtained  from  Fig.  (2a),  is  9.7  jisec  as  compared  to 
the  theoretical  valu?  of  10.42  n sec.  The  small  difference  is  attributed  to  the  fact  that 
r0  and  T\  depend  on  the  volume  of  the  numerical  integration.  It  should  be  noticed  that 
in  the  numerical  calculation  the  integration  has  been  carried  inside  a  cylinder  of  radius 
300  cm  and  of  half-height  200  cm,  while  in  the  theoretical  work  the  integration  has  been 
carried  out  over  the  whole  space.  To  obtain  some  insight  into  this  effect,  let  us  look  at 
the  diffusion  through  a  cylindrical  conducting  shell  of  infinite  length,  where  the  region  of 
the  numerical  integration  is  limited  within  a  perfect  conductor  located  at  a  radius  rj .  The 
diffusion  time  is  equal  to  rp  [l  —  (a/ri)2],  where  a  is  the  radius  of  the  conductor  and  rp 
is  given  by  Eq.  (21).  It  is  apparent  that  the  diffusion  time  decreases  as  the  volume  of  the 
integration  region  decreases,  and,  therefore,  the  inductance  of  the  system  decreases.  As 
another  example  consider  a  toroidal  conducting  shell,  with  a  region  of  integration  confined 
within  a  perfect  conductor  located  at  77  =  r)i  and  therefore  the  radial  positions  r  are 
limited  to  the  values  between  r\  mtn  =  6tanh(r7i/2)  and  rj  moi  =  bcoth(r)i/2).  All  the 
theoretical  results  obtained  above  remain  the  same  except  for  the  L/R  time  in  Eq.  (20a) 


11 


that  is  replaced  by  2 T& 


<n!T2  _  2  _  (V^oY 

\  i  max  ) 


Again,  the  L/R  time  (but  not  the 


,2ri  max  ' 

times  Ti,  r2)  decreases  as  the  volume  of  the  integration  region  decreases  and  so  does  r^. 
In  the  example  given  above  the  region  for  r  <  rimtn  is  excluded,  while  in  the  numerical 
results  this  region  is  included.  For  this  reason,  the  comparison  of  the  numerical  results 
is  made  with  the  theoretical  results  for  the  whole  space  rather  than  a  finite  region  of 
integration.  The  largest  discrepancy  between  theory  and  numerical  computation  is  in  the 
field  index  for  small  times.  This  could  be  due  to  the  large  time  step  (0.5  nsec)  that  was 
used  in  the  numerical  integration  or  to  the  fact  that  the  field  index  is  computed  from  the 
second  derivative  of  the  stream  function  rAg  and  this  function  is  not  accurately  computed 
initially  by  the  code.  The  same  argument  is  applicable  to  the  smaller  discrepancy,  for 
short  times,  in  the  quantity  $o/{nroBzo).  Over  all,  however,  the  agreement  is  fairly  good 
between  theory  and  numerical  computation  with  the  TRIDIF  code. 


III.  Experimental  Results 

This  section  briefly  summarizes  the  measurements  on  the  time  delay  of  the  vertical 
magnetic  field  caused  by  the  vacuum  chamber  of  the  device  and  compares  the  experimental 
results  with  the  predictions  of  the  theory  and  that  of  the  computer  code. 

The  NRL  modified  betatron  comprises  three  different  external  magnetic  fields;  the 
betatron  field  that  is  a  function  of  time  and  is  responsible  for  the  acceleration  of  the 
electrons,  the  toroidal  magnetic  field  that  varies  only  slightly  during  the  acceleration  of 
the  electron  ring  and  the  strong  focusing  field  that  also  has  a  very  weak  time  dependence. 
These  three  fields  have  been  described  previously.  However,  for  completeness  we  briefly 
describe  in  this  paper  the  vertical  field  and  the  vacuum  chamber. 

The  betatron  magnetic  field  controls  mainly  the  major  radius  of  the  gyrating  electron 
ring  and  is  produced  by  a  set  of  18  air  core  (see  Fig.  3),  circular  coils  connected  in  series. 
Sixteen  main  coils  and  four  trimmers  are  used  to  generate  a  field  configuration  with  a  flux 
condition  equal  to  1.92.  Their  total  inductance  is  approximately  343  /jH.  The  coils  are 
powered  by  an  8.64  mF  capacitor  bank  that  can  be  charged  up  to  17  kV.  At  full  charge, 
the  bank  delivers  to  the  coils  a  peak  current  of  about  45  kA.  The  current  flowing  through 
the  coils  produces  a  sinusoidally  varying  field  having  a  quarter  period  of  2.6  msec  and  its 
amplitude  on  the  minor  axis  at  peak  charging  voltage  is  2.1  kG.  Immediately  sifter  the  peak 
the  field  is  crowbarred  with  a  4.5  msec  decay  time.  The  temporal  profile  of  the  vertical 
field  is  shown  in  Fig.  4. 
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The  flux  condition  and  field  index  are  adjusted  by  two  sets  of  trimmers  that  are 
connected  in  parallel  to  the  main  coils.  The  current  through  the  trimmers  is  adjusted  with 
series  inductors.  Typically  10  -  15%  of  the  total  current  flows  through  the  trimmers. 

The  100  cm  major  radius,  15.2  cm  inside  minor  radius  vacuum  chamber  has  been 
constructed  using  epoxy  reinforced  carbon  fibers.  The  desired  conductivity  is  obtained  by 
embedding  a  phosphor  bronze  screen  inside  the  body  of  the  graphite  as  shown  in  Fig.  5. 
The  graphite  is  2.5  mm  thick  and  has  a  surface  resitivity  of  26.6  mfl  on  a  square.  The 
screen  has  250  x  250  wires  per  inch  and  is  made  of  40  ^m  diameter  wire  with  an  equivalent 
surface  resistivity  of  12.8  mfl  on  a  square.  The  calculated  resistance  for  the  entire  vacuum 
chamber  is  57  mfl.  The  measured  D.C.  resistance  of  the  toroidal  vacuum  vessel  is  68  ±  2 
mfl.  The  outside  surface  of  the  chamber  is  covered  with  a  6.3  mm  thick,  epoxy  reinforced 
fiberglas  layer.  Figure  6  is  a  photograph  of  the  vacuum  chamber. 

This  novel  construction  technique  has  several  attractive  features,  including  control¬ 
lable  resistivity  and  thus  magnetic  field  penetration  time,  high  stiffness  and  tensile  strength, 
high  radiation  resistance  (up  to  500  Mrad)  and  low  outgassing  rate  (~  10-8  Torr/sec-cm2). 

The  vertical  magnetic  field  is  monitored  with  a  small  magnetic  probe  that  is  located  on 
the  minor  axis.  The  probe  has  been  wound  of  gauge  35  copper  wire,  it  has  approximately 
600  turns,  its  measured  inductance  is  6.2  mH  and  its  measured  internal  resistance  is  43  fl. 
The  output  of  the  probe  is  fed  to  a  passive  integrator  with  a  time  constant  of  20.5  msec. 

To  improve  the  time  response  of  the  probe  the  50Q  terminator  at  the  input  of  the 
integrator  has  been  omitted.  Oscillations  in  the  output  signal  have  been  avoided  by  locating 
the  oscilloscope  approximately  2  m  away  from  the  probe,  inside  an  aluminum  housing 
cladded  with  mu  metal. 

Results  from  the  measurements  are  shown  in  Fig.  7.  The  lower  trace  is  the  output 
of  the  probe  when  it  is  located  inside  the  vacuum  chamber  and  the  upper  trace  is  the 
output  of  the  probe  after  the  chamber  is  removed  (vacuum  field).  It  is  apparent  from  the 
oscillogram  that  after  approximately  60  nsec  the  two  signals  are  parallel  and  the  lower  is 
delayed  from  the  upper  by  34  n sec. 

For  the  parameters  of  the  vacuum  chamber,  and  when  the  flux  condition  ^o/tttqBzo 
is  equal  to  1.92,  i.e.,  for  t  =  0.04  both  theory  and  the  computer  code  predict  a  time  delay 
of  37  jxsec.  The  less  than  10%  difference  between  the  two  delay  times  is  probably  within 
the  uncertainty  of  the  measurement. 

In  addition,  we  broke  the  electrical  continuity  of  the  vacuum  chamber  by  unbolting 
the  joint  of  two  adjacent  sectors.  The  probe  was  placed  on  the  minor  axis  with  its  axis  of 
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symmetry  on  the  vertical  plane  that  passes  through  the  symmetry  plane  of  the  gap.  The 
intermediate  trace  in  Fig.  7  shows  the  output  of  the  probe  for  this  measurement. 

It  is  expected  that  at  the  symmetry  plane  of  the  gap  the  magnetic  field  to  be  identical 
to  the  vacuum  field.  However,  due  to  various  constraints,  the  diameter  of  the  probe  and  the 
gap  width  are  approximately  equal  and  thus  the  measured  field  is  lower  than  the  vacuum 
field. 

IV.  Conclusion 

The  diffusion  of  the  magnetic  field  through  a  toroidal  conducting  shell  has  been  studied 
under  the  assumption  of  a  small  aspect  ratio.  The  external  magnetic  field  can  have  an 
arbitrary  field  index  and  magnetic  flux  on  the  minor  axis  of  the  torus.  The  diffused  field, 
field  index,  magnetic  flux  and  wall  current  were  computed  analytically  and  compared 
with  the  numerical  results  from  the  TRIDIF  code.  Three  time  constants  determine  the 
evolution  in  time  of  the  diffusion  process,  namely,  the  L/R  time  to,  the  diffusion  time  Tj, 
and  T2  =  ti/2.  The  delay  time  depends  linearly  on  To,  Ti  and  also  on  the  flux  condition  of 
the  external  field.  The  analytic  delay  time  was  larger  to  that  computed  from  the  TRIDIF 
code.  The  difference  is  attributed  to  the  finite  volume  of  the  region  of  integration  in 
the  numerical  computation  which  causes  the  inductance  of  the  system  to  be  smaller.  In 
general,  the  agreement  between  the  theoretical  and  numerical  results  is  quite  good. 

A  measurement  of  the  delay  time  in  the  toroidal  chamber  of  the  NRL  modified  be¬ 
tatron  gave  a  delay  time,  approximately  equal  to  34  fisec,  i.e.,  less  than  10%  smaller  from 
the  theoretical  value  of  37  /zsec. 

In  conclusion,  we  have  seen  that  the  diffusion  process  in  the  toroidal  conducting  shell 
is  much  more  complicated  than  that  in  a  conducting  cylinder.  Therefore,  the  results  for  a 
cylinder  cannot  be  generalized  to  apply  to  a  toroidal  device.  The  case  of  the  torus  should 
be  investigated  on  its  own  merit. 
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Table  I. 


Parameters  for  the  diffusion  theory. 

Torus  major  radius,  ro  (cm)  100 

Torus  minor  radius,  a  (cm)  16 

Chamber  wall  thickness,  dw  (cm)  2 

Torus  wall  conductivity,  a  (mho/cm)  14.666 

Amplitude  of  sinusoidal 

external  magnetic  field,  bo  (Gauss)  1267 

Period  of  external  magnetic  field,  T  (//sec)  104 

Field  index  of  external  magnetic  field,  n  0.53 

Flux  condition  parameter 

of  external  magnetic  field,  e  -0.05 
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Toroidal  shell  and  coordinate  systems  used  in  the  analysis 


VERTICAL  MAGN  FIELD  (GAUSS) 
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Fig.  2  —  (continued)  Temporal  profiles  from  analysis  of  the  (solid  line)  and  computer  code  (dashed  line), 
(a)  vertical  (betatron)  field,  (b)  field  index,  (c)  flux  condition,  and  (d)  wall  current. 
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Fig.  3  —  Elevation  of  the  experiment 
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Fig.  4  —  Vertical  (betatron)  field  vs.  time 
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VACUUM  CHAMBER  WALL  COMPOSITION 


See  Detail  A 


30°  Sector  of  Chamber 

100cm  Major  Radius 
15  cm  Minor  Radius 
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(250  X  250-7^—,  1.6  milwire) 
Inch 
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Low  Outgass  Epoxy 

Atmosphere 


Detail  A 

Wall  Construction 

Fie  5  —  30°  sector  of  the  vacuum  chamber  and  details  of  wall  construction 


fie  6  Photograph  of  the  epoxy  reinforeed-carbon  fiber  vacuum  e.iamber 


Hi:  "  Measured  vertical  magnetic  field  vv  nine  of  the  minor  axis  (r  =  100  cm.,  z  =  0)  of  the 

torus  Bon. mi  trace:  Inside  the  complete  vacuum  chamber  with  a  gap.  Top  trace:  After  the 
removal  ot  the  chamber  (vacuum  Held). 
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